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ABSTRACT 

In massive primordial galaxies, the gas may directly collapse and form a single central massive object 
if cooling is suppressed. H2 line cooling can be suppressed in the presence of a strong soft-ultraviolet 
radiation field, but the role played by other cooling mechanisms is less clear. In optically thin gas, 
Lyman-a cooling can be very effective, maintaining the gas temperature below 10'' K over many orders 
of magnitude in density. However, the large neutral hydrogen column densities present in primordial 
galaxies render them highly optically thick to Lyman-a photons. In this paper, we examine in detail 
the effects of the trapping of these Lyman-a photons on the thermal and chemical evolution of the 
gas. We show that despite the high optical depth in the Lyman series lines, cooling is not strongly 
suppressed, and proceeds via other atomic hydrogen transitions. At densities larger than ~ 10^ cm~^, 
collisional dissociation of molecular hydrogen becomes the dominant cooling process and decreases 
the gas temperature to about 5000 K. The gas temperature evolves with density as T cx p^"'*"', with 
7eff = 0.97 — 0.98. The evolution is thus very close to isothermal, and so fragmentation is possible, 
but unlikely to occur during the initial collapse. However, after the formation of a massive central 
object, we expect that later-infalling, higher angular momentum material will form an accretion disk 
that may be unstable to fragmentation, which may give rise to star formation with a top-heavy IMF. 
Subject headings: atomic processes - molecular processes - stars: Population HI - cosmology: theory 
- dark ages, reionization, first stars 



1. INTRODUCTION 

Supermassive black holes with masses M > 10*- 
10^ M© are known to have e xisted at very earl y times in 
the history of the universe (iJiang et al.| [200f). prompt- 
ing one to ask how black holes of such a size could have 
formed so quickly. They may have started life as stellar- 
mass objects, the compact remnants of the first genera- 
tion of stars, and they may have grown to their present 
size by accretion. However, forming supermassive black 
holes in this way presents a number of difficulties. The 
first stars are thought to have had masses of order lOOM© 
or mo re, making them strong sources of ionizing radi- 
atioii (lAbel et al.l 120021: iBromm fc LarsonI 120041: iGloveii 
[20051 : lYoshida et al.l I2008n . The dark matter halos in 
which these stars formed were small, with masses of only 
lO^M©, and models of the effects of the ionizing radiation 
produced by a primordial star within one of these small 
halos show that it readily photo-evaporates the bulk of 
the gas from the halo, leaving, little gas available to 
be accreted by the black hole (poh nson & Bromm 20Q3; 
IMilosavlievic et al.|[200llAlvare^'et al.ll2009l) . Although 
accretion is likely more efficient in the 5 — 6 cr peaks 
that harbor the observed supermassive black holes, it is 
stil l constrained b y the Eddington limit. As discussed 
bv IShapird ([2005I) . seed masses of at least 10^ Mq are 
required for a n MHD disk or a standard thin disk. Simi- 
larly, fSchlcich er et al.l pOlO.) found that high Eddington 
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ratios and unusually low mass-energy conversion efficien- 
cies are needed for stellar progenitors. 

Previous studies have therefore examined the pos- 
sibility of forming these large seed black holes 
by d i rect gravitational collap se (lEisenstein &: Loebl 
'1995; KoushiaoDas eLaU 12004 iBegelman et al.l 120061 : 
.Spaans fc Silk 200^. Gas falling into a halo with a 
mass M > 5 X 10^[(1 -I- z)/10]-3/2 Mq (corresponding 
to a virial temperature Tvir > lO** K) will be shock- 
heated up to Tvir, and so will become hot enough to 
cool via the electronic emission lines of atomic hydro- 
gen. The fate of the gas in one of these so-called 'atomic 
cooling' halos will then depend on its subsequent ther- 
mal evolution. If the gas can cool efficiently as it col- 
lapses, significantly lowering both its temperature and 
its Jeans mass, then it is very difficult to prevent it 
from fragmenting, a nd hence forming stars rather than a 
mass i ve black hole (IBromm fc Loebll2003l : lOmukai et al.l 
[2008t iClark et al.l 120091). On the other hand, if cool- 
ing remains inefficient and the effective equation of state 
for the gas is stiffe r than is otherm al, then fr agmenta- 
tion is suppressed (iLi et al.l ^2005; Lodato fc N ataraianl 
|2006( ) and the formation of a massive black hole is a 
plausible outcome. In gas that has been enriched with 
heavy elements and dust from the first generation of 
stars, it app ears that cooli ng and fragmentation cannot 
be avoided (jOmukai et al.l [2008) . On the other hand, 
in primordial gas the only effective low-temperature 
coolant is molecular hydrogen (H2), and if enough of 
this can be destroy ed by photodissociation, then cooling 
can be suppressed IB egelman et al." 2006l: iDiikstra et al.l 
[20081 IRegan fc Had mclt 2009; Shang,eialjl2009D. The 
soft-UV background needs to be strong enough to sup- 
press H2 formation until densities of ~ lO'"^ cm~'^ are 
reached, following which collisional dissociation will sup- 
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press the H2 abundance for gas temperatures greater 
than a few 1000 K. The required background field de- 
pends on the chemical model, the spectral shape and 
the potential presence of additional electrons created 
in sho ck front s. It is of t he order J21 = 10^ — 10 ^ 
([Omik ai 2001: iBromm fc Lo eb 2003: iShang "eral1[2009l ) . 
These values are much larger than estimates of the mean 
strength of the backg round, which range from J21 ^ 0. 1 
pohnson et all 120081) to J21 ~ 40 piikstra et al.ll2Q08D . 
However, recent work has shown that the radiation back- 
ground is highly inhomogeneous, owing to the strong 
clus tering of the first generation of star-fo rming galax- 
ies ()Diikstra et all 120081: lAhn et alj [20090 . In regions 
close to these gala xies, the ultrav i olet b ackground can 
be very large, and iDiikstra et al.l (|2008[ ) have recently 
demonstrated that enough massive primordial galaxies 
with a local value of J21 > 10'^ may exist. In addition, 
early reionization in local patches of the universe may 
suppress star f ormation in so me primordial halos in a 
similar fashion ()Johnsonll2009fl . They may thus remain 
pristine, but grow in mass until they can collapse. 

In the absence of H2 cooling, emission from the Ly- 
man series lines of atomic hydrogen, especially Lyman-a, 
plays a central role in regulating the temperature of the 
collapsing gas. Hydrodynamical models of the collaps- 
ing gas typically assume that the Lyman series lines re- 
main effectively optically thin, allowing the optically thin 
form of the cooling rate to be used. These models find 
that the temperature of the collapsing gas decreases with 
increasing density, albeit very slowly, typically reaching 
T - 6000-7000 K at densities of order n ~ 10^ cm'^. 
However, the large hydrogen column densities present in 
these protogalaxies produce very large optical depths in 
the Lyman series lines, and so the assumption that the 
optically thin form of t he cooling rate can saf ely be used 
is highly questionable. ISpaans &: Silkl ()2OO60 examined 
the effects of Lyman-a photon trapping within an ana- 
lytic framework for the evolution of the gas and showed 
that in the absence of additional cooling mechanisms, it 
gives rise to a temperature evolution described by 

dlogT 0.5- 7/18 -Bn^/is 

~ rfk^ ^ log Cni/2 + SnVis ' 

where p is the mass density, T the temperature, n 
the number density, B ^ 0.5 — 0.1 cm^/^ and C 
10-36^^/(10'^ Mq) cm3/2, .^^-^j^ ^.^^ ^Yie halo mass. 
For halo masses in the range 10^-10^ M©, this yields 
7eff — 1 0.01 — 0.5, indicating an effective equation of 
state that is stiffer than isothermal, and hence a temper- 
ature that increases with increasing density. This result 
suggests that not only is the neglect of the effects of 
photon trapping a poor approximation, but it may also 
lead to qualitatively incorrect results, as it is crucial for 
fragmentation whether the effective equation of state is 
harder than isothermal. 

In this paper, we re-examine the issue of Lyman-a pho- 
ton trapping, using a considerably more detailed treat- 
ment than in Spaans & Silk (2006). In particular, we ad- 
dress the issue of whether there are other cooling mech- 
anisms that can compensate for Lyman-a if the latter is 
strongly suppressed. 

2. METHOD 



To follow the thermal and chemical evolution of the 
collapsing gas, we make use of a one-zone treatment in 
which the form of the density evolution is prescribed in 
advance. One-zone models are widely used for study- 
ing the chemistry and thermodynamics of primordial 
or low metallicity gas (e. g . lOiriukai 2001: Omukai"et al.l 
20051: IC azaux fc Soaansl 120091: IGl" over fc SavinI I200a 
Schleicher et al.ll2009D , as they enable one to model the 
chemistry in great detail and to include a large number 
of different cooling processes, without the computational 
efficiency concerns inherent to a three-dimensional treat- 
ment. 

We extend the one- zone model develope d by 
IGlover fc Sa^ ((20091) and lSchleicher et al.l (l2009l) to in- 
clude effects that become important at higher temper- 
atures and in the presence of Lyman-a trapping. Near 
T ^ 10^ K, H" form ation cooling m ay contribute to the 
overall cooling rate (|Omukaill200lD . We assume that a 
typical electron undergoing radiative attachment to form 
H~ has an energy of order /c^T, where is kB Boltzmann's 
constant, and hence write the H~ formation cooling rate 
as 

A-H- = kn-nnnckBT, (2) 

where fcjj- is the H" formation rate coefficient, nn is 
the atomic hydrogen number density, and rio is the elec- 
tron number density. To model the effects of Lyman-a 
trapping, we include different level populations as sep- 
arate species in the code. In our model, we consider 
energy levels up to n = 5. For the first excited state, 
we distinguish between the 2s and the 2p states, as only 
decays from the 2p state to the ground state will produce 
Lyman-a photons; two-photon decays from the 2s state 
will produce continuum photons that will not be trapped. 
The other states are considered as averages over the an- 
gular momentum quantum numbers. Radiative decay 
rates Aij, collisional excitation and de-excitation rates 
Cij, collisional ionization and three-body recombination 
rates, and collisional transition rates between the 2s and 
2p state are all adopted from Omukai (2001). We define 
the transition rate Rij from level i to level j as 

(1 + Q,y) + C„-, i>j (3) 

^ij ~ -^jil^cscJiQ ji ~t- Cij ^ i <^ J, (4) 

9i 

where /3osc,ij is the escape probability for the i ^ j tran- 
sition, gi — denotes the statistical weight of level i 
and Qij — Jcont,^ /(2/iffj), where Vij is the frequency 
of the transition and Jcont.ij the average intensity of the 
background radiation field at this frequency. The spec- 
tral shape of the backgrou nd radiation fi eld is taken to 
be that of a 10^ K (IGlover fc SavinI I2009D . and the nor- 
malization is set by fixing the value of J21. We assume 
that UV photons more energetic than 13.6 eV are ab- 
sorbed by atomic hydrogen in the intergalactic medium, 
and we do not consider the effects of an extragalactic X- 
ray background. For radiative or collisional transitions 
to the first excited state, we assume that the reaction 
products will be distributed according to the statistical 
weights of the 2s and 2p states. However, our results are 
not sensitive to this assumption, as collisional transitions 
between these two states occur rapidly in the conditions 
of interest, and so the ratio between the 2s and 2p level 
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populations is always very close to equilibrium. 

For the escape probability Pcsc,ij, we adopt the expres- 
sion 



1 - exp(-Ti, 



-exp(-^tphAcoii) 



(5) 



where Tij is the optical depth at line centre of the i ^ j 
transition, tcoii the collapse time of the gas, and tph is 
the photon diffusion time, i.e. the time required for a 
photon to diffuse out of the optically thick gas, and /3 a 
geometrical factor. Even small deviations from spherical 
symme try lead to escape along a preferred direc tion and 
/3 = 3 (jPiikstra et all [20061 : ISi^ aans fc Silkl[2006h . In our 
model, the collapse time is the free-fall time, corrected 
by a factor t hat takes into accou i it the effective equatio n 
of state (see lOmukai et all 120051: iSchleicher et al.ll2009[ ). 
For most line transitions, tph <C tcoU, such that the ex- 
ponential factor is negligible. For these lines, we follow 
lOmukail (pOOO) and assume that the dominant contri- 
bution to the optical depth comes from material within 
one local Jeans length, and that the density, temperature 
etc. do not vary significantly on this scale. This is justi- 
fied provided that the collapse is close to isothermal, and 
that the initial mass of the collapsing gas is comparable 
to the Jeans mass, as in this case the density and veloc- 
ity profiles of the gas will c ome to resemble the Larson - 
Penston similarity solution ( jLarsoniri969l : lPenstonlll969D . 
which has just these properties. For Lyman-a pho- 
tons and other direct transitions to the ground state, 
we compute the optical depth r and the pho ton dif- 
fusion timescale following ISpaans fc Silkl (|2006( ) , in or- 
der to take into account the significant line broaden- 
ing that occurs during the diffusion process. For the 
Lyman-a line, this yields a timescale iph — L{aT2iY^^ /c 
with the natural-to-thermal line-width a and an optical 
depth T21 = 1.04 X lO'^^NuT^^-^ , with L the distance 
to the edge of the halo, and TVh the column density of 
atomic hydrogen. For the Lyman-a line, the natural- 



to-thermal line-width is taken as a = 4.7 x 10 ^T. 



-1/2 



([Spaans &: Silk|[2006[) , while we correct this value for the 
increased lifetime in case of other line transitions to the 
ground state. To calculate the column density, we as- 
sume that the density profile of the gas scales with ra- 
dial distance r as r~^' ^, as indicated by p revious numer- 
ical simulations fe.g. iWise fc Abeill2007D . The cooling 
functions for the hydrogen lines are then evaluated based 
on the escape probability, the level populations and the 
strength of the background radiation. 

3. RESULTS AND CONCLUSIONS 

The temperature evolution that we find with our ap- 
proach for several different values of J21 is plotted as a 
function of density in Fig. [TJ Even for J21 = 1, the tem- 
perature remains much higher during the collapse than 
in the radiation-free case, and never drops significantly 
below 1000 K. 

As previously noted, much higher values of J21 can be 
obtained locally, in the presence of a luminous neighbou r 
within a distance of 10 kpc or less (iDijkstra et al.ll2008l ). 
As we increase J21, we find systematically higher temper- 
atures at each density, and for J21 = 100, we find that 
the effects of H2 collisional dissociation start to become 
important. In this case, II2 cooling is marginally effective 
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Figure 1. Temperature evolution as a function of density for dif- 
ferent values of J21. The thin dotted lines indicate lines of constant 
Jeans mass. 
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Figure 2. Fractional abundances for electrons, H2 and the hydro- 
gen level populations, plotted as a function of density, for the case 
where J21 = 1000. The values of the hydrogen level populations 
depend on the assumed size of the protogalaxy; values are plotted 
for virial radii of 500 pc (thin lines) and 3 kpc (thick lines). 
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Figure 3. The main cooling functions for J21 = 1000 as a function 
of density. Lines are black if the contribution is not affected by 
Lyman a trapping. For the hydrogen lines, blue lines correspond 
to a system with virial radius 0.5 kpc (~ lO'^ Mq), and red lines 
to a virial radius of 3 kpc (~ 2 X 10^ Mq). The H2 collisional 
dissociation cooling rate is corrected for the effect of H2 heating. 

at densities n < 10^ cm~'^, but at higher densities, colli- 
sional dissociation destroys the II2, causing a sharp rise 
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in the temperature. Such a radiation field may be present 
in a fraction of 10~^ — 10^^ of all atomic cooling halos 
(jPiikstra et all [20081 ) . For J21 = W and J21 = 10'*, we 
find, in common with previous studies, that H2 cooling 
never becomes important, and the gas temperature re- 
m ains high. Nevertheles s, in contrast to the predictions 
of iSpaans fc Silkl ()2006[ ). we find that the temperature 
does decline with increasing density, although it never 
drops below 5000 K over the range of densities examined 
here. With 7off — 0.97 — 0.98, we find that the effective 
equation of state of the gas is just slightly softer than 
isothermal. We note, though, that this estimate provides 
just a lower limit for the value of 7eff, which is strictly 
defined as d\ogT/d\ogp at constant entropy. 

In order to understand why the gas is still able to cool, 
we explore the case with J21 = 1000 in more detail in 
Fig. [2] and Fig. [3l In Figure [21 we show the evolu- 
tion of the fractional abundances of free electrons and 
H2 molecules, as well as the fractional level populations 
of the 2s, 2p and 3 states of atomic hydrogen (denoted 
henceforth as H(2s), H(2p) and H(3)). In Figure [3l we 
show the contributions that the main cooling processes 
make to the total cooling rate. The Lyman-a optical 
depth depends on the size of the protogalaxy, and we 
consider two examples, one with a virial radius of 500 pc 
(corresponding to a halo mass of ^ 10^ Af©), and a sec- 
ond with a virial radius of 3 kpc (corresponding to a halo 
mass of 2 X 10^ Mq). In each case, the abundance of 
the 2p state is initially very small, but increases roughly 
linearly with increasing density, owing to the increasing 
importance of Lyman-a photon trapping. However, col- 
lisional excitation to the 11(3) state eventually comes to 
dominate over radiative decay from the 2p state, follow- 
ing which the abundance of the 2p state remains roughly 
constant. 

The abundances of the H(2s), H(2p) and H(3) states 
are clearly larger in the case with the larger virial ra- 
dius, as Lyman-a photons are trapped more effectively, 
directly boosting the population of the 2p state. The 
population of the 2s state is then increased due to colli- 
sional coupling. The same holds true for the H(3) state. 
However, the electron and H2 abundances are not af- 
fected. Despite these changes in the hydrogen level pop- 
ulations, the total cooling rate does not change signifi- 
cantly. Cooling due to the 2p-ls line decreases, but this 
is balanced by additional cooling from the 2s-ls and 3- 
2 transitions, due to the increased populations of these 
states. At densities higher than 10^ cm~'^, collisional 
dissociation of H2 becomes the dominant cooling process 
(even after correcting for H2 formation heating) and cools 
the gas to about 10^ K. H~ formation cooling is found 
to contribute significantly in a broad range of densities. 

Although we find that the effective equation of state of 
the gas may be softer than isothermal, it is nevertheless 
true that 7cff ~ 1 throughout the collapse. Moreover, the 
temperature evolution of the gas is surprisingly similar 
to that obtained in previous three-dimensional simula- 
tions that assume optically thin Lyman-a cooling. The 
outcome of the initial collapse may thus be similar to 
that found in previous simulations, i.e. the formation of 
a single massive bound object at the ce ntre of the halo 
(see e.g. [W ise et al. 2008; "Begclman & Shlosmanl [2001 
iRegan fc Haehneltl 120091: IShang et al.l (20091). However, 



it is also highly likely that gas that falls in later, with 
higher angular momentum, will begin to build up an 
accretion disk surrounding this object. The high gas 
temperatures imply a rapid accretion flow onto such a 
disk, which may quickly become gravitationally unsta- 
ble. If the gas in the disk is also able to cool effectively, 
which our results suggest will be the case even if Lyman- 
a cooling is completely discount ed, then we would expect 
it to fragment and form stars (|Levin|[2007[ ). The latter 
may be drawn to the center by d ynamical fric tion and 
merg e with the central black hole (jPevecchi fc Volonterl 
|2009( ). Owing to the high temperatures, we would ex- 
pect the stars formed in the disk to accrete more rapidly 
than standard Population HI stars, by a factor of ten or 
more, and hence their final masses may be very large, 
M 100 — 500 A/q. However, recent works also indicate 
the possibility of quasi-stable cold self-gravitating accre- 
tion disks that do not fragment, which c ould feed the 
central object during the further ev olution ([Levine et al.l 
[20081 IBegelman fc Shlosim^ [20091) . 
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